Vortex solitons in photonic crystal fibers 
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We demonstrate the existence of vortex soliton solutions in photonic crystal fibers. We analyze the 
role played by the photonic crystal fiber defect in the generation of optical vortices. An analytical 
prediction for the angular dependence of the amplitude and phase of the vortex solution based on 
group theory is also provided. Furthermore, all the analysis is performed in the non-paraxial regime. 



Spatial soliton solutions in ID periodic structures 
have been theoretically predicted ,1] and experimentally 
demonstrated in nonlinear waveguide arrays 1^. More 
recently, similar structures have been observed in opti- 
cally induced gratings 0, 0] and, for the first time, in 
2D photonic lattices 1^. Modeling has been tradition- 
ally performed using the discrete nonlinear Schrodinger 
equation (NLSE)[llg, valid only in the so-called tight- 
binding approximation. However, the accurate study of 
nonlinear solutions in optically-induced lattices requires 
the resolution of the NLSE with a continuous periodic 
potential model |^ Q |^ . Discrete optical vortices have 
been also predicted in the discrete NLSE and in con- 
tinuous models jS] and, recently, experimentally observed 
in optically- induced 2D lattices 0. 

A different scenario appears in a 2D photonic crystal 
in the presence of a defect. An example of such a system 
is a photonic crystal fiber (PCF). A PCF is a silica fiber 
that possesses a regular array of air-holes extending along 
the entire fiber length (see Fig^^a)). It is an axially 
invariant 2D photonic crystal with a central defect (the 
PCF core, where guidance occurs). It has been recently 
proven that a PCF acting as a 2D nonlinear photonic 
crystal can support and stabilize fundamental solitons 
[iflj. Unlike fundarnental solitons embedded in perfectly 



periodic potentials a- distinguishing feature of these 

systems is that no power gap is needed to generate the 
nonlinear localized solution. In this letter we go a step 
further and give numerical evidence of the existence of 
vortex solitons in PCF's. 

When one considers that silica can have a nonlinear 
response, electromagnetic propagation in a PCF is given 
by the following nonparaxial equation (V • E w 0): 



k[) = ujq/c the vacuum wavenumber. The linear refrac- 
tive index profile function nQ(x,y) is 1 in the air-holes 
and equals risiiica in silica, whereas the nonlinear index 
profile function n2{x,y) is different from zero only in 

silica {n\,a,,,) 



3X(smca)/(2^ocno(siiica)))- We search 



for monochromatic (or quasi-monochromatic) stationary 
electric field solutions with well-defined constant polar- 
ization: E(a;,y, z,t) « u0(x, y)e'('^^"""*l That is, we 
analyze the following equation: 



L{m = i3^^ Lm = L^ + LNLm), (2) 

where Lq = ^1 + klnl{x,y) and Lnl{\(I)\) = 
k'^n\{x,y)\4>\'^ stand for the linear and nonlinear parts 
of the differential operator L(|0|), respectively. Due to 
the geometry of the air-holes distribution, both no{x,y) 
and n2{x,y) are invariant under the action of the Ce^ 
point-symmetry group. This group is constituted by dis- 
crete 7r/3-rotations {Rtt/^) plus specular reflections with 

respect to the x and y axes: y ^ 



[V^ + klinlix) + n2(f)|E|2)]E = -d^-E/dz^ 



(1) 



being the 2D-transverse Laplacian operator and 



-y and x 

equivalently, 9 ^ —9 and 6 — > tt — 9, in polar coordinates. 

Our aim is to look for nonlinear solutions of Eq.jpil be- 
yond the fundamental solitons reported in Ref . ilfll . Our 
search strategy will be based on group theory arguments, 
which will be enlightening in the classification and un- 
derstanding of new solutions according to the sixth-fold 
symmetry of the system. In the linear case, this approach 
has been adopted to study the degeneracy of solutions in 
a PCF Here, we will use it to obtain analytical 

expressions for the angular dependency of solutions. Ac- 
cording to group theory, since Lg is invariant under the 
Cqv group, all its eigenfunctions have to lie on discrete 
representations of this group 12] . Solutions of the vortex 
type belong to one of the two two-dimensional represen- 
tation of Cqv. a pair of fimctions {<j>i,4>*) {I = 1,2) be- 
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FIG. 1; (a) Schematical representation of a PCF. (b)-(d) 
Amplitudes of several vortices for increasing values of 7 in 
a PCF with A = 31/xm and a = 6/xm (A = 1064 nm): (b) 
7 = 0, 45 X 10"^ (c) 7 = 0, 95 X 10"^ (d) 7 = 1, 75 x 10"^ 



longing to these representations have the following trans- 



formation properties: {(t>h(t>i) i^'' 4>U f-*^ (t>i) e = e^^/'^, 

^ 0/) and ^ These 

symmetry constraints impose the form of these functions: 



0f(r,0)exp[z0f(r,0)] 1 = 1,2, 



(3) 



where (r, 0) is a scalar function, characterized by 
(^"(r,0 + 7r/3) = (l)''{r,d) and 0"(r, -6*) = (f>^{r,7T-0) = 
(p^{r,6), and (jf{r,9) is a pseudoescalar function char- 
acterized by (j)P{r,e + 7r/3) = (j)P{r,e) and ^^''{r, -6*) 
0P(r,7r - 6*) = -(^P(r,6l). Both 0" and (jf are peri- 
odic functions of 0, so they can be expanded in Fourier 
series in cos(6n6') and sin(67i0) {n £ N). Reflection 
symmetry forces (/)*((/)^) to depend on the cosine (sine) 
terms only: i^f (r, 6) = C'ini'f) cos(6n0) and </if (r, 9) = 
Y^n bin{r) sin{&nO). 

In the nonlinear case, functions 4>i of the form given 
by Eq.Q can satisfy self-consistency, according to 
group theory. The full operator L[\(j)i\) = Lo{nQ) + 
LNL{r''(l)'^;n2) is invariant under the Cqv group because 
Lo, r' , (t)i and 712 are all Cqv invariants. Thus L{\(j)i\), like 
in the linear case, provides the representation where (jji 
lies on. However, group self-consistency is not sufficient 
to ensure that the resolution of the nonlinear equation 
^ with the ansatz Q provides a nontrivial solution in 
all cases, since Eq.Q always admits the cj) = solu- 
tion. We solve Eq. (|2Jl by means of the Fourier iterative 
method previously used in Ref. 10] to find soliton solu- 
tions in PCF's. These solutions also fulfill the group self- 



consistency condition and they belong to the fundamen- 
tal representation of Cg^. The important difference now 
is that we restrict ourselves to a different representation 
space; we search for nonlinear solutions in the Z = 1,2 
representations of Cqv The solution and its conjugate 
0;* represent a vortex and an anti- vortex soliton of order 
I, respectively. 
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FIG. 2: Effective index of a family of vortex solutions as a 
function of 7 (solid line) . Same for a family of fundamental 
solitons (dashed line), as in Ref..lO!]. The shadow region corre- 
sponds to the conduction band, constituted by Bloch modes, 
of the 2D photonic cladding with A = 31 fim and a = 10 fim 
(A = 1064 nm). 

We have simulated large-scale PCF's (in order to pre- 
vent silica breakdown |lO|) with a lattice period, or 
pitch, A = 31 fim, and several air-hole radius at a fixed 
wavelength A = 1064 nm (a convenient quasi-continuum 
source). We define a dimensionless nonlinear param- 



eter as 7 



2(silica) 



P/Aq (Aq is an area parameter: 



Aq = 7r(A/2)^; P is the total power carried by the so- 
lution) and perform calculations for increasing values of 
7. We discover that a family of optical vortex solutions of 
the type given by the ansatz ||2J) is found. The amplitudes 
of several solutions (with / = 1) for increasing values of 
7 are represented in Fig^b)-(d). According to group 
theory (Eq.|(21)), these amplitudes have to be scalar func- 
tions, thus showing full invariance under C^y, and they 
have to vanish at the origin, as seen in Figmb)-(d). Be- 
sides, they become gradually narrower as the nonlinearity 
increases. It is interesting to plot the effective index of 
these solutions (rieff = P/ko) versus 7, as shown in FigI21 
The value of riefi increases as 7 increases, accordingly to 
the narrowing of solutions depicted in Fig^b)-(d). For 
comparison, we have also included the neff(7) curve of 
the fundamental soliton family for the same PCF struc- 
ture, as in Ref. [loj . Both lie on the upper forbidden band 
of the perfectly periodic cladding structure. Here it can 
be clearly envisaged the role of the central defect. Notice 
that there is no threshold 7 (power) to generate a nonlin- 
ear vortex soliton, i.e, there is a continuum of solutions in 
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7 starting from the linear (TE or TM) mode. Unlike in 
perfectly periodic structures, the defect can eliminate the 
presence of a threshold power If we consider struc- 
tures where the linear mode is not present |l3| , a thresh- 
old 7 is necessary to generate the nonlinear solution in 
the forbidden band. By playing with the geometric pa- 
rameters of the PCF, it is possible to tune this threshold 
value. In all cases, this value is much lesser than in per- 
fectly periodic structures, which can be important for 
experimental purposes. 

Even more interesting is the phase behavior of these so- 
lutions. We show in Fig|3|a typical phase profile of a PCF 
vortex calculated at a fixed radius. It shows strict accor- 
dance with the group theory prediction for the phase: 
arg(0O + (/.f (r, 6) = W + Y.n=i ^Ur) sin(6n0). The 
phase of these vortex solutions differs from that of an or- 
dinary vortex by the presence of the pseudo-scalar func- 
tion extra term: the group phase. Besides the existence 
of the standard linear behavior in 0, the group phase 
provides an additional sinusoidal dependence on the an- 
gle (with period determined by the group order: 7r/3 for 
order 6) not present in ordinary vortices. 




FIG. 3: Phase of a vortex with I — 1 
represent both the total phase a.Tg{(j>i) 
group phase arg(<^;) — W (dashed line). 



Sit r = 21 /im. We 
(solid line) and the 



In order to check the stability of a vortex soliton solu- 
tion we have to consider z-dependent perturbations. 
This implies solving the non-paraxial equation for 
the perturbed field (f)' = + Scj). In terms of group 
theory, {(f)i,(j)*) constitute a basis of the I representa- 
tion of Cqv, so we can consider two types of perturba- 
tions around a vortex solution (/);: diagonal perturba- 
tions, in which (j)' remains in the one-dimensional sub- 
space spawned by 0;, and non-diagonal, where the per- 
turbation takes (f)' out of this subspace. An initial con- 
dition of the form (f>' — re** exp[— (r/rg)^] is an exam- 
ple of diagonal perturbation for a vortex {I = 1). Non- 
paraxial evolution shows that a vortex is stable under 
diagonal perturbations. In Fig^J we show the stabiliza- 

,1/2 



asymptotic vortex state. Diagonal perturbations gener- 
alize the concept of radial perturbations to the discrete 
symmetry case. We have also simulated non-diagonal 
perturbations. In this case, an oscillatory instability oc- 
curs, as seen in Fig0f inset). Non-diagonal perturbations 
correspond to azimuthal ones in the radially symmetric 
case. Despite the presence of azimuthal instabilities, no 
collapse is detected due to the non-paraxial nature of 
evolution 
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FIG. 4: Axial evolution of {nef;{z)} for a diagonal perturba- 
tion. Inset: two snapshots of the asymptotic behavior of a 
non-diagonal oscillation instability (notice the different posi- 
tions of the peak intensity). 

Summarizing, we propose a novel way of generating 
vortex solitons, based on PCF's, in which the presence of 
the defect plays a crucial role. Our analytical approach, 
based on group theory, is completely general and applies 
to any system owning a discrete invariance, no matter is 
periodic or non-periodic, linear or nonlinear. Finally, the 
use of a non-paraxial equation permits to analyze new 
phenomena beyond the paraxial approximation. 
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